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The count rate of Weakly Interacting Massive Particle (WIMP) dark matter can- 
didates in direct detection experiments experiences an annual modulation due to the 
Earth's motion around the Sun. In the standard isothermal halo model, the signal 
peaks near June 2nd at high recoil energies; however, the signal experiences a phase 
reversal and peaks in December at low energy recoils. We show that this phase reversal 
may be used to determine the WIMP mass. If an annual modulation were observed 
with the usual phase (i.e., peaking on June 2nd) in the lowest accessible energy recoil 
bins of the DAMA, CDMS-II, CRESST-II, EDELWEISS-II, GENIUS-TF, ZEPLIN-II, 
XENON, or ZEPLIN-IV detectors, one could immediately place upper bounds on the 
WIMP mass of [103, 48, 6, 97, 10, 52, 29, 29] GeV, respectively. In addition, detec- 
tors with adequate energy resolution and sufficiently low recoil energy thresholds may 
determine the crossover recoil energy at which the phase reverses, thereby obtaining 
an independent measurement of the WIMP mass. We study the capabilities of various 
detectors, and find that CRESST-II, ZEPLIN-II, and GENIUS-TF should be able to 
observe the phase reversal in a few years of runtime, and can thus determine the mass 
of the WIMP if it is 0(100 GeV). Xenon based detectors with 1000 kg (XENON and 
ZEPLIN-IV) and with energy recoil thresholds of a few keV require 25 kg-yr exposure, 
which will be readily attained in upcoming experiments. 

1 Introduction 

Extensive gravitational evidence suggests that a dominant fraction of the matter in our 
Galaxy is nonluminous, or dark. Although the identity of this dark matter is currently 
unknown, it may consist of Weakly Interacting Massive Particles (WIMPs). Numerous 
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experiments aimed at the direct detection of these WIMPs are currently being developed. 
These experiments generically measure the nuclear recoil energy deposited in a detector 
when an incident WIMP interacts with a nucleus in the detector. An important signature of 
halo dark matter in direct-detection experiments is the annual modulation induced by the 
Earth's motion with respect to the halo jT]. In the standard model of the dark halo, in which 
the velocity distribution of the WIMPs is a Maxwellian distribution truncated at the escape 
speed of the Galaxy, the modulation of the WIMP interaction rate at high nuclear recoil 
energies is in phase with the motion of the Earth with respect to the halo, peaking when 
the relative speed is maximal (June) and reaching a minimum when the relative motion 
is minimal (December). For low energy recoils, however, it is well known that that the 
interaction rate experiences a phase reversal, peaking in December [21 El El E] • This phase 
reversal occurs below a particular crossover recoil energy, Qc- In this paper we examine 
what information about the mass of the WIMP can be obtained from this phase reversal. 
We emphasize that this phase reversal could constitute an important signature of a WIMP 
flux, and we also point out that because the crossover recoil energy is a function of the WIMP 
mass, the observation of the phase of annual modulation immediately places an upper limit 
on the allowed WIMP mass. 

Furthermore, if WIMP direct detection experiments with sufficient energy resolution 
determine the crossover energy Qc, the WIMP mass could be independently measured using 
only observations of the phase of the annual modulation. 

This paper is organized as follows. In section 2, we briefly review the basics of WIMP 
direct detection experiments and the isothermal model of the dark halo. In section 3, we 
explore the dependence of the annual modulation on energy recoil, and revisit the phase 
reversal that occurs at low recoil energies. Finally, we demonstrate how this observation 
may be exploited to place limits on the WIMP mass and to obtain estimates of the mass 
directly. We emphasize that observation of this phase reversal will constitute an unambiguous 
signature of an extraterrestrial WIMP flux. 

2 WIMP Direct Detection Experiments 

More than twenty collaborations worldwide are presently developing detectors designed to 
search for WIMPs. Although the experiments employ a variety of different methods, the 
basic idea underlying WIMP direct detection is straightforward: the experiments seek to 
measure the energy deposited when a WIMP interacts with a nucleus in the detector ^T] . If 
a WIMP of mass scatters elastically from a nucleus of mass ttiat, it will deposit a recoil 
energy oi Q = (m^f^/mAr)(l — cos6'), where 

rrir = m^rriN / {m^ + mj^) (1) 



2 



is the reduced mass, v is the speed of the WIMP relative to the nucleus, and 9 is the 
scattering angle in the center of mass frame. We may compute (following, e.g., |T71 ITHl IT^) 
the differential detection rate, per unit detector mass {i.e., counts/day/kg detector/keV 
recoil energy) associated with this process, 

';^ = ^F\Q)T{Q,t) (2) 

where ph = 0.3 GeV/cm'^/c^ is the halo WIMP density, ctq is the total nucleus- WIMP inter- 
action cross section and F{Q) is the nuclear form factor for the WIMP-nucleus interaction 
that describes how the effective cross section varies with WIMP-nucleus energy transfer. We 
consider here only spin independent interactions, wherein the target nucleus can be approx- 
imated as a sphere of uniform density smoothed by a gaussian ^3], and the resulting form 
factor is, 

j7(g) _ 3[sin(Qri) -gricos(gri)] ^„q2,2/2 

where ri = (r^ — bs'^Y^'^, s ^ 1 fm, and, 

r ~ [0.91(M/GeV)^/^ + 0.3] x lO^^^cm (4) 

is the radius of the nucleus. For spin-dependent interactions, the form factor is somewhat 
different but again -F(O) = 1. In general, the form factor must be evaluated for each detector 
nuclei. For a more extensive discussion, see fn\ [T| IT^ ITHl IT7j. 
For purely scalar interactions, 

477? 

aO,scalar = -^[Z fp + {A - Z) f^f . (5) 
TT 

Here Z is the number of protons, A — Z is the number of neutrons, and fp and /„ are the 
WIMP couplings to nucleons. For purely spin-dependent interactions, 

c^o,spin = (32/7r)G^/i2AV( J + 1) . (6) 

Here J is the total angular momentum of the nucleus and A is determined by the expectation 
value of the spin content of the nucleus (see [TT[ [TJ [13 CHI El)- 

For the estimates necessary in this paper, we take the WIMP-nucleon cross section^ 
CTp = 7.2 X 10""^^ cm^, and take the total WIMP-nucleus cross section to be ^ 



^This is the cross section with nucleons at zero momentum transfer as discussed in Eq. (7.36) of |17j . 
^In most instances, /„ ~ fp so that the fohowing equation results from Eq.©. 
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where the rrirp is the proton- WIMP reduced mass, and A is the atomic mass of the target 
nucleus. 

Information about halo structure is encoded into the quantity T{Q,t), 

T{Q,t)=r^dv, (8) 

where fd{v) is the distribution of WIMP speeds relative to the detector, and where v^i^i = 
{QrriN /2mlY^'^ represents the minimum velocity that can result in a recoil energy Q. To 
determine T{Q,t), we must have a model for the velocity structure of the halo. 

2.1 The isothermal halo 

In what follows we will consider dark matter detection experiments without directional ca- 
pabilities; such experiments are sensitive only to the WIMP flux integrated over the entire 
sky. The most frequently employed background velocity distribution is that of a simple 
isothermal sphere. In such a model, the galactic WIMP speeds with respect to the halo obey 
a Maxwellian distribution with a velocity dispersion ah, 

/ 3 \3/2 / \ 

= 4, «^exp^-^j (9) 

where v is the WIMP velocity relative to the Galactic halo, and where we have performed an 
angular integration over the entire sky. We take the velocity dispersion of our local halo to 
be ah = 270 km/s. This is the local velocity distribution; it is expected to vary with spatial 
position throughout the Galaxy. 

As pointed out by Drukier, Freese, and Spergel P, and studied (in the context of an 
isothermal halo) by Freese, Frieman, and Gould Earth-bound observers of the dark halo 
will see a different, time dependent velocity distribution as a result of the relative motion of 
the Earth with respect to the Galaxy. To take this into account, we estimate the velocity of 
the Earth with respect to the Galactic halo. Neglecting the ellipticity of the Earth's orbit 
and the non-uniform motion of the Sun in right ascension, we may write the speed of the 
Earth with respect to the the dark WIMP halo as 

Veh = \l + Vlf^ + 2VesVsh COs(27rt - (ph) (10) 

where the phase (ph = 2.61 ± 0.02 corresponds to June 2nd ± 1.3 days (relative to = 
on January 1). The speed of the Earth with respect to the Sun is Ves = 29.8 km/s, and 
the speed of the Sun with respect to the halo is Ves = 233 km/s. More precise expressions 
for the motion of the Earth and Sun may be found in |UI2]. Translating the distribution in 
Eq. iniinto the distribution as seen by an earthbound detector, fd{veh), and integrating the 
resulting distribution in Eq. |Sl one obtains 
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With Th{Q, t), one can now compute the recoil energy spectrum associated with the isother- 
mal halo, as well as the annual modulation of that spectrum that results from the motion 
of the Earth around the sun. For the case of an isothermal halo, one generally expects that, 
above a critical recoil energy, the annual modulation will be in phase with the motion of 
the Earth with respect to the halo, peaking near June 2nd. In the following section, we will 
demonstrate the dependence of the phase of the modulation on the observed recoil energy. 

3 The Phase Reversal at Low Recoil Energies 

The annual modulation of the WIMP interaction rate |2], and in particular the phase 
reversal at low recoil energies [3 13 El , has been studied in detail by several authors. Here 
Fig. ^plots the expected differential event rate for a WIMP of mass 65 GeV at recoil energies 
of 35 keV and 10 keV. We note that there is a 180 degree phase shift between the two signals. 
The origin of this effect may be understood by considering a first order approximation to 
Eq. ^2 We note that the velocity of the Earth with respect to the WIMP halo, Veh, given 
in Eq. ITUlmav be written, for Ves <^ Veh, as 



Veh ~ Vh + ^hit). 



where we define, 



Vk = ivi + vt,)y\ (12) 
= 7I^X;XW^c°s(27rt-0;,). (13) 

Using this approximation. We may expand this expression about e(t) = 0, an approximation 
that is well justified because eh(t)/rih ^ 1 for all t. With the following definitions, 

= erff^%^^V (14) 



2' 



= exp(- ^^—; ''"^ , (15) 



we find that, to first order in e. 



5 




Figure 1: The annual modulation of the WIMP differential detection rate for a WIMP with 
a mass of 65 GeV, in a ^'^Ge detector, for recoil energies of 35 and 10 keV. Note the 180 
degree phase shift between the two signals. 
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8 Et + El 



c^hVh 



(16) 



We see that the differential detection rate will be modulated by e(t) oc cos(27rt — (ph), as 
expected. We are particularly interested in the behavior of the amplitude of the oscillation. 
If we rewrite Th{Q, t) of Eq. [TBI as 



T(Q, t) = B + A cos(27rcjt - (j)h) 
we identify the amplitude of oscillation as 



;i7) 



A{Q) 



VesVsh 




8 E^ + Ef^ - Xf^ 



O-hVh 



For A{Q) > 0, the modulation of the WIMP differential event rate is in phase with the 
motion of the Earth with respect to the halo, and will peak at June 2nd as usual. For 
A{Q) < 0, however, the modulation of the event rate moves 180 degrees out of phase with 
terrestrial motion, peaking in December rather than June. Values of A{Q) are plotted in 
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Fig. Elfor < Q < 100 keV, for a WIMP mass of 30, 50 and 70 GeV in a ^^Ge detector. 
In general, for a WIMP mass and nuclear mass rriN there will exist a nonzero, crossover 
recoil energy, Qc, at which A{Qc) = 0, and below which A{Q < Q^) < 0. That is, for 
Q < Qc the WIMP interaction rate will peak in December rather than June. We emphasize 
that these approximations, while illustrative, give rise to errors of a few percent; in practice, 
Eq. [TT] should be used for all computations. 




10 20 30 40 50 60 70 80 90 100 
Recoil Energy Q (keV) 



Figure 2: The amplitude of modulation A{Q) as a function of recoil energy Q for WIMP 
masses of 30, 50, and 70 GeV in a "^^Ge based detector. For each WIMP mass there exists 
a particular recoil energy Qc, at which the first order annual modulation of the WIMP 
differential event rate vanishes. For Q < Qc, we find that A[Q < Qc) < 0, and the differential 
WIMP interaction rate peaks in December, rather than June. 

The crossover recoil energy, Qc is a function of the mass of the WIMP, as well as the mass 
of the target nuclei in the detector. We can estimate the critical recoil energy, below which 
we observe the phase reversal. Fig. Olplots the crossover recoil energy as a function of WIMP 
mass, for Nal, ''^Ge, ^^^Xe, CaW04 and AI2O3 based WIMP direct detection experiments. 
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Collaboration 


Material 


Qthresh (keV) 


Md (kg) 


CDMS-II 




10^ 


5 


CRESST-II 


CaW04 


~ 1 


10 


DAMA 


Nal 


22 


100 


EDELWEISS-II 


73Ge 


20 


38 


GENIUS-TF 


^^Ge 


1 


40 


XENON 


LXe 


4 


1000 


ZEPLIN-II 


LXe 


10 


40 


ZEPLIN-IV 


LXc 


4 


1000 



Table 1: Properties of various dark matter detection experiments including: the detector 
material, the detector energy recoil threshold, Qthresh, and the total detector mass, M^). In 
the case of future experiments, the quantities are projected. 

3.1 Placing an Upper Bound on the WIMP Mass 

The observation of an annual modulation signal can be used to constrain the WIMP mass. 
For example, if a direct-detection experiment observes an annually modulated signal peaking 
in June at an energy recoil bin Qo, we may place an upper limit on the WIMP mass because 
we know that the crossover energy Qc niust have been less than Qq, or a phase reversal would 
have been observed. Because the crossover energy is a monotonically increasing function of 
the WIMP mass, this places an upper limit on the WIMP mass. 

We will describe the upper bounds one can obtain in a variety of dark matter detectors. 
Of course, the experiments must have sufficient exposure to see the annual modulation in 
order to obtain these bounds. 

As a concrete example, consider the DAMA experiment, which reports an annual modu- 
lation in the event rate in its Nal detector in the 22-66 keV recoil energy bins [211. Although 
other experiments rule out much of this region of parameter space j2Z] , if one were to believe 
the DAMA results, the fact that annual modulation in DAMA peaks in June in this energy 
bin implies, from Fig. El that the maximum WIMP mass consistent with the DAMA data 
is m < 103 GeV. 

This same logic may be applied to other dark matter experimental collaborations that 
may be sensitive to the phase reversal. In particular, in this paper we will consider the 
various versions of the CRESST CDMS 26j, EDELWEISS (27!, XENON GENIUS 
[5U] . and ZEPLIN jSI] experiments. Table I collects the detector target materials, mass, 
and recoil energy thresholds associated with these experiments. In Table II, we present the 
maximum WIMP mass consistent with an observation of a June-peaking annual modulation 
in the lowest accessible recoil energy bins of the respective experiments, as determined using 
the arguments in the preceding paragraphs. 



8 



Collaboration 


Maximum allowed 
WIMP mass (GeV) 


CDMS-II 


48 


CRESST-II 


6 


DAMA 


103 


EDELWEISS-II 


97 


GENIUS-TF 


10 


XENON 


29 


ZEPLIN-II 


52 


ZEPLIN-IV 


29 



Table 2: The maximum WIMP masses consistent with observing a June-peaking annual 
modulation signal in the lowest accessible energy recoil bins of the respective dark matter 
experiments. 

3.2 Detecting the Phase Reversal 

For WIMP direct detection experiments with sufficient energy resolution and sufficient ex- 
posure, it may be possible to determine the crossover recoil energy itself, thereby making 
an independent measurement of the WIMP mass. If one determines the recoil energy at 
which the modulation phase changes sign. Fig. El may be used to estimate the WIMP mass. 

We now explore whether or not the phase reversal itself could be practically observed 
in the context of current and near-future direct detection experiments. In order to extract 
a small annual signal from the background, we require a large signal to noise ratio so that 
the fluctuations are small compared to the desired signal. This problem, in the context of 
WIMP direct detection experiments, has been discussed by Hasenbalg 0. Because the count 
rate for direct detection experiments is so low, extended exposure times will be required to 
reliably detect a small modulation. 

We write the total signal as a function of time as. 



S{t) = / —dQ (19) 
= So{Qi,Qf) + SUQi,Qf)cos{ut) + 0{Sl) (20) 

where 5*^ is the amplitude of modulation, 5*0 is the unmodulated count rate, and dR/dQ is 
defined as in Eq. |21 These quantities will depend on the limits of integration, Qi and Qf, 
and may be be written in terms of the differential event rate evaluated at its maximum in 
June, Sj, and the rate at its minimum, in December, Sd, 
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Figure 3: The crossover recoil energy Qc as a function of WIMP mass for Nal, ^^Ge, 
^^^Xe, CaW04, and AI2O3 detectors, assuming a halo velocity dispersion of 270 km/s. The 
observation of an annual modulation signal peaking in June at a given recoil energy neces- 
sarily places an upper limit on the allowed WIMP mass. For example, the observation of a 
June-peaking annual modulation in the 22 keV energy recoil bin implies an upper limit on 
the WIMP mass of about 103 GeV for a Nal based experiment. 



Sm = \[Sj-SD\ So = \[Sj + Sd\. (21) 

If Sra ^ •S'o, the theoretical signal to noise ratio may be written as, 

is/n) = -^^^^lElLVMf. (22) 
/SoiQ.,Qf) 



where M is the total detector mass and T is the total exposure time. As a reasonable 
criterion for distinguishing the modulation signal from the noise, we require that the s be at 
least 2a greater than the statistical uncertainty. This amounts to requiring (s/n) = 2. For 
low count rates, large MT will be required to achieve this minimum signal to noise ratio. 

In order to observe a phase reversal in a particular WIMP detector with a given WIMP 
mass, an annual modulation peaking in December must be detected in the energy range 
Qthresh to Qc, and a modulation peaking in June must be observed in the energy range Qc to 
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Qmax, where Qthresh IS the recoil energy threshold of the detector, and Qmax is the maximum 
recoil energy accessible to the detector. In Fig. HI we plot the minimum exposure (MT)min 
(in kg-yrs) necessary to observe a phase reversal for several current and future experiments, 
requiring a minimum signal to noise ratio {s/n) = 2. We should note that we are assuming 
perfect energy detector resolution, and that these exposure estimates therefore represent 
lower limits to the required exposure in actual experiments. 
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Figure 4: The minimum exposure (in kg-yrs) required to observe the phase reversal in the 
WIMP-nucleus differential event rate at the 2a level, as a function of WIMP mass for various 
detectors. 

We can determine the viablity of observing the phase reversal and hence obtaining an 
estimate of the WIMP mass in various detectors. We note that, using this phase reversal 
technique, the experiments are most sensitive to WIMP masses between 80 GeV to a few 
hundred GeV, which is the most interesting mass range. Accelerators (LEP-200) place lower 
bounds on the neutralino mass of roughly 45 GeV. The upper limit for sensible supersym- 
metric dark matter particles is roughly 1.5 TeV. Recent upper bounds on the Constrained 
Minimal Supersymmetric Model are tighter |18|IT^ I^. Hence the most sensitive mass range 
in the detectors covers a compelling slice of parameter space. 
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3.3 Summary of Results 

As emphasized above, Fig. IH presents the total exposure time required to observe the phase 
reversal as a function of the WIMP mass, for various detectors. 

The CRESST-II experiment, which employs calcium tungstate (CaW04) crystals, should 
be able to detect the phase reversal. In particular, for WIMP masses near 80 GeV, only two 
years of total exposure would be required to observe the phase reversal. 

For the ^^Ge based detectors of CDMS and CDMS-II, the phase reversal is not likely to 
be observed. CDMS-II has projected exposures of roughly 2,500 kg-days of data [2H]- For 
the most favorable WIMP masses, nearly 30,000 kg-days of data are required to observe the 
phase reversal using germanium detectors with a recoil energy threshold of 10 keV. 

The EDELWEISS-II experiment, which also employs ^^Ge, the relatively high energy 
recoil threshold of the experiment renders it unlikely to observe the phase reversal. At best, 
for WIMP masses near 400 GeV, the experiment requires ~ 340 kg-yrs of exposure. With 
its current mass, this will necessitate nearly 9 years of continuous running time. 

The GENIUS-TF experiment, a test facility for the projected GENIUS experiment, is 
currently the best germanium-based platform for observing the phase reversal [20] • Requir- 
ing only 10 kg-yrs of exposure for WIMP masses ~ 80 — 120, the experiment could observe 
the phase reversal with as little as a year of runtime. The germanium based experiments 
have the advantage over other experiments in that they are sensitive to a larger range of 
WIMP masses due to the relative flatness of the {MT)^[^ curves in Fig. Other germanium 
based experiments not reflected in the plot, such as the 0.2 kg HDMS experiment [32] do 
not have sufficient exposure to observe the reversal. 

The DAMA experiment, which is not represented in Fig. |31 utilizes low mass Nal de- 
tectors with a high energy recoil threshold of 20 keV. The DAMA collaboration will not be 
able to see the phase reversal, despite its already extensive, 60,000 kg-day exposure. For 
the most optimistic WIMP masses, DAMA would require nearly 1000 kg-yrs of exposure to 
observe the WIMP phase reversal. 

Liquid xenon experiments that are currently underway, such as the 40 kg ZEPLIN-II, will 
be sensitive to the phase reversal [3T]. With a recoil energy threshold of 10 keV, the ZEPLIN- 
II experiment will require only 40 kg-yrs of exposure for WIMP masses near 150 GeV. For 
these WIMP masses, it may observe the phase reversal with only a year of continuous running 
time. 

The upcoming, large scale LXe detectors such as ZEPLIN-IV and XENON with masses 
on the order of 1000 kg and energy recoil thresholds of a few keV are in a excellent position 
to observe the phase reversal. LXe detectors with energy recoil thresholds on the order of a 
few keV require 25 kg-yrs exposure for WIMP masses in the range 100-150 GeV. Hence, such 
experiments could observe a phase reversal with as little as a year of continuous runtime. 

As seen in Fig. HI a disadvantage of xenon based detectors is that the minimum required 
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exposure is low only for a narrow span of WIMP masses. Outside of the 100-200 GeV 
mass range, the required exposure rapidly becomes larger. The other detectors, particularly 
CRESST-11 and GENIUS-TF, have flatter curves, i.e., the minimum required exposure does 
not change rapidly as a function of WIMP mass. However, future 1 tonne xenon detectors, 
such as the XENON and ZEPLIN-IV experiments, will have the advantage of sheer size, so 
that they may reach the required exposures on reasonable timescales. 

4 Conclusions 

The reversal of phase at low recoil energy constitutes an unambiguous signature of an inci- 
dent WIMP flux. The usual annual modulation, while an important indicator of a WIMP 
halo, could conceivably be mimicked by Earth bound physics (such as annually modulated 
temperature fluctuations in the detector). It is difficult, however, to imagine non-WIMP 
scenarios that could engender a phase reversal only at low energy recoils. In addition, the 
critical recoil energy depends in a well understood way on the WIMP mass, so that an 
independent estimation of the WIMP mass can be obtained by locating the critical recoil 
energy, Qc, below which the phase reversal is observed. Furthermore, because the crossover 
recoil energy is a monotonically increasing function of the WIMP mass, the observation of an 
annual modulation in the WIMP signal will necessarily imply an upper hmit to the WIMP 
mass. 

If a June-peaking annual modulation were observed in the lowest accessible energy recoil 
bins of the DAMA, CDMS-II, CRESST-II, EDELWEISS-II, GENIUS-TF, ZEPLIN-II or 
XENON/ZEPLIN-IV detectors, one could immediately place upper bounds on the WIMP 
mass of [103, 48, 6, 97, 10, 52, 29] GeV, respectively, provided the experiments have sufficient 
exposure to observe the annual modulation at all. 

The observation of a phase reversal would provide independent confirmation of the WIMP 
mass. The CaW04 based detectors of CRESST-II have the lowest required exposure, on the 
order of 10 kg-yrs, and would require only 2 years of running time to observe the phase 
reversal. With 40 kg of detector, GENIUS-TF and ZEPLIN-II could observe the phase 
reversal with as little as a year of runtime. We note that the 1000 kg Xenon based detectors, 
ZEPLIN-IV and XENON, operating with recoil energy thresholds on the order of a few keV 
provide an excellent means of detecting the phase reversal in the differential rate for WIMP 
masses on the order of 100 GeV, requiring only 1 year of running time. Hence, should the 
annual modulation be observed, its phase reversal will provide a means to obtain the WIMP 
mass in these detectors. 

The analysis of this paper has assumed an isothermal halo model. Galaxies, however, are 
thought to form hierarchically, with small substructures condensing first and subsequently 
merging into progressively larger structures. It is therefore reasonable to think that the dark 
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halo of our galaxy may be populated with clumps or streams of dark matter. Indeed, work by 
Stiff, Widrow, and Frieman suggests that there is a high (C(l)) probability that residual 
substructure exists at the solar radius, with a contribution that is an additional few percent 
of the local Halo density. In fact, as discussed by ISHESl, the Sagittarius dwarf galaxy, which 
is being tidally disrupted by the Milky Way, is pouring dark matter down upon the solar 
neighborhood, at the level of an additional (0.3-25)% of the local Halo density. The peak 
of the annual modulation and value of the crossover recoil energy will differ in halo models 
containing substantial clumps or streams of dark matter, or characterized by anisotropic 
velocity distributions jSH]. In these alternate halo models, the phase and amplitude of the 
annual modulation are generically altered, and the impact on the crossover recoil energy 
must therefore be computed on a case by case basis. Such a study will be the subject of 
future efforts. 

We thank J. Frieman, H. Nelson, P. diStefano and L. Stodolsky for helpful remarks. We 
acknowledge the DOE and the Michigan Center for Theoretical Physics at the University of 
Michigan for support. 
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